Enhanced decoherence in the vicinity of a phase transition 
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We study the decoherence of a spin- 1/2 induced by an environment which is on the verge of a 
continuous phase transition. We consider spin environments described by the ferromagnetic and 
antiferromagnetic Heisenberg models on a square lattice. As is weU known, these two dimensional 
systems undergo a continuous phase transition at zero temperature, where, the spins order sponta- 
neously. For weak coupling of the central spin to these baths, we find that as one approaches the 
transition temperature, critical fluctuations make the central spin decohere faster. Furthermore, 
the decoherence is maximal at zero temperature as signalled by the divergence of the Markovian 
decoherence rate. 
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The dynamics of a real system is determined not only 
by its internal Hamiltonian but also by its environment. 
Due to the coupling to the numerous environmental de- 
grees of freedom, an initial pure quantum state of the 
system evolves into an incoherent mixture. This pro- 
cess, decoherence, is a major hurdle in the construction 
of quantum computers where sufficiently long coherence 
times of the qubits (the basic units of quantum infor- 
mation) are fundamental requisites. The recent spate of 
experimental work on qubits has generated a great deal 
of interest in the question of decoherence induced by dif- 
ferent environments. The most studied environments are 
bosonic baths modeled as baths of harmonic oscillatorsl. 
More recently, due to the numerous experimental realiza- 
tions of qubits which involve real spin- 1/2 objects^, there 
has been an increasing focus on spin baths as a primary 
source of decoherence. 

The rich and varied physics of spin systems make spin 
baths fundamentally interesting. Spin baths comprising 
non-interacting spins have been among the first studied 
models to understand the decoherence process-. How- 
ever, since the physical origin of the decoherence of a 
system lies in the dynamical fluctuations of the bath de- 
grees of freedom to which it couples, we expect the re- 
sulting decoherence to reflect the non-trivial nature of 
the fluctuations induced by the interactions in the bath. 
From this perspective, the vicinity of a continuous phase 
transition is especially interesting given the existence of 
critical fluctuations. Clearly, these divergent fluctuations 
are expected to have dramatic consequences for the de- 
coherence. Though some authors have studied spin bath 
models which exhibit continuous phase transitions, the 
effect of critical fluctuations has been occulted either be- 
cause of the purely mean field nature of the models^ or 
due to the temperature regimes considered^i^. 

In this Letter, we address the issue of the decoherence 
of a central spin weakly coupled to a spin bath which is on 
the verge of a phase transition. We model the spin bath 
as a two-dimensional system of Heisenberg spins with 
nearest neighbour interactions on a square lattice. We 
consider both ferromagnetic and antiferromagnetic inter- 



actions. In both cases, the spin bath undergoes sponta- 
neous symmetry breaking at Tc = to a magnetic phase. 
These spin bath models are interesting for our purpose 
as their critical fluctuations are well understood^. Stud- 
ies on fluctuation free mean fleld models indicate that 
the decoherence time scale increases as temperature is 
lowered and the system undergoes a transition'*. More- 
over, this result is in accord with the conventional wis- 
dom that decoherence is minimal at T = and increases 
with increasing thermal fluctuations. In what follows, we 
will show that contrary to the above scenario, as one ap- 
proaches the critical point at zero temperature, critical 
fluctuations result in a faster decoherence of the central 
spin with the decoherence being maximal at zero tem- 
perature. 

We consider a central spin-1/2 coupled isotropically via 
a hyperfine like contact interaction to a spin-5 bath with 
uniform nearest neighbor interactions. For simplicity, we 
assume that the central spin has no internal dynamics. 
The total Hamiltonian describing the central spin and its 
environment is given by 



H = Hi + Hb 

= cr • ^ AjS'i - • 5j 

i (ij) 
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where cr are the Pauli matrices, S{ the spin operators 
of the bath spin at site i and (...) denotes summation 
over nearest neighbor spin pairs of a square lattice with 
periodic boundary conditions. For the case of an elec- 
tronic spin interacting with nuclear spins through the 
hyperfine contact interaction, the coupling constant A; 
is simply related to the electron envelope wave function 
at the site i. Antiferromagnetic interactions (J < 0) 
and ferromagnetic interactions (J > 0) lead to Neel or- 
der and ferromagnetic order respectively, at T = 0. To 
study the decoherence for weak coupling to the spin bath 
and a finite temperature T, we use the resolvent operator 
approach^'^'^*' in conjunction with the Schwinger boson 
technique. To simplify the calculation, we assume a fac- 
torizable initial density matrix for the composite system. 
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= p{Q)®pB where p(0) and pB oc exp{—HB/T) denote, 
respectively, the central spin state and the thermal equi- 
librium of the bath. We use units ks ~ h = 1 throughout 
this paper. 

Since the total Hamiltonian ([T]) is rotationally invariant 
and the finite temperature phase of the bath is isotropic, 
the central spin Bloch vector at time t is related to its 
initial value as given by 



(2) 



where r is a scalar function of the time t and (...) denotes 
the average over the density matrix of the composite sys- 
tem. Due to the coupling to the bath, the function r 
vanishes in the long time limit. To determine r{t), it 
is useful to write the reduced density matrix p{t) of the 
central spin as 
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where is a real positive number, Tr^ denotes the par- 
tial trace over the bath degrees of freedom and L is the 
Liouvillian corresponding to the total Hamiltonian i?, 
i.e., CA — [H,A\ for any operator A. Using the decom- 
position of the density matrix p{t) in the basis of 2 x 2 
matrices {I,aa} and the projection operator technique 
explained in Refs. and [TO, we obtain from ([3]) the ex- 
pression ^ where 
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The self-energy S is given by 



= iTr [a^Trs [Ci{z - Q C Q)-^ CipB<Ta\ 



(4) 
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where ctq, is any Pauli matrix and Ci is the Liouvillian 
corresponding to the interaction Hamiltonian Hi. The 
projection operator Q is defined by its action on any op- 
erator A as QA = A — TTBiA)pB- We note that for 
any arbitrary Hamiltonian, the dccohcrence of the cen- 
tral spin cannot be described by a single self-energy, the 
full time evolution of the state p{t) is given by a 4 x 4 
matrix function of the complex variable z. From this 
4x4 matrix, it can be shown that the asymptotic be- 
havior of p(t) is in general characterised by two different 
times, a decoherence time and a relaxation time which 
determines the time for thermal equilibration of the cen- 
tral spinii. For the Hamiltonian ([1]), the time-scales for 
relaxation and decoherence are the same. However, for 
the generic case of the central spin having its own inter- 
nal dynamics, these two times are in principle different. 
This generic case is studied in the last part of the paper. 

The decoherence in the weak coupling regime is deter- 
mined by the lowest order contribution of the interaction 
Hamiltonian Hj to the self-energy S. The first term of 
this expansion is given by (O with C replaced by the bath 
Liouvillian Cb- Using the properties of the Pauli matri- 
ces, this second-order contribution to the self-energy can 



be written as 

S2(^) = -^TrE^i^J / dte'^'ReiSiit) ■ S^)b (6) 
^ ij 

where Si{t) = exp{iHBt)Sieyip{—iHBt) and {A)b = 
Ti {pb A) denotes the thermal average of any bath op- 
erator A. Neglecting higher order contributions to S in 
is equivalent to the Born approximation^^. It can be 
shown Ifi that in the weak coupling limit, the expression 
(HI can be simplified to obtain 



Inr(t) ~ 

TT 



,E^2^r,iE; (7) 



where r2(i;) = -ImT.2{E+iO+). Eq.jT]) describes the de- 
coherence at all time scales. It yields the usual quadratic 
decrease r{t) ~ 1 — J dE T2 {E)/2tt for very short times 
and Markovian decay r{t) ~ exp[— r2(0)i] for asymptotic 
times. The full time evolution of the central spin is de- 
termined by the dynamic spin correlation function of the 
bath via the relations ([7]) and ([6]). Typically, dynamical 
correlations are rather difficult to calculate for spin sys- 
tems. In the following, we use the successful Schwinger 
boson mean field theory to evaluate the dynamical spin 
correlation of the bath described by Hb- We will show 
that the critical fluctuations in the bath lead to a diver- 
gence of r2(0) at the transition, implying a faster than 
exponential asymptotic decay of r{t). This divergence is 
merely a reflection of the divergence of the underlying 
correlation time in the bath at the phase transition. 

In the Schwinger representatiorti^, at every site i, the 
spin operators S\ are replaced by two bosonic opera- 
tors aj and 6; with the constraint a\a-^ + b\b.^ — 2S. 
In the ferromagnetic case, the correspondence relations 



read 2S? 



b]b: and 5*;"^ = alb: . In the antifer- 



romagnetic case, due to the possibility of Neel ordering 
of spins, two sub-lattices must be distinguished. On one 
of the sub-lattices, the spin and boson operators are re- 
lated as above, while on the other sub-lattice the cor- 
respondence relations take the form 25'? = — a|aj 

and ~ both ferromagnetic and antiferro- 

magnetic cases, the Hamiltonian ^ can be interpreted 
as describing a spin-1/2 coupled to a boson bath. The 
resulting problem is nonetheless very different from the 
standard spin-boson model for the following reasons: i) 
the central spin couples to two species of bosons ii) the 
bath bosons interact with each other and are subject to 
constraints which conserve the number of bosons. In the 
Schwinger boson mean field theory, the local constraints 
on the bosons are replaced by a global constraint via a 
uniform chemical potential and the boson Hamiltonian 
derived from Hb is studied using a Hartree-Fock mean 
field schemei^iii^'i^. We now present analytical results for 
the decoherence obtained within this theory in the low 
temperature regime T <C \J\S. 

For the ferromagnetic Heisenberg model (J > 0), the 
mean-field approximation for Hb yields, up to a constant 
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energy, = I]k^k(akak + ^k\) where ajj and by, 

are the Fourier transforms of a; and 6;, respectively. The 
magnon dispersion is given by 



Wk = 2JQ(2 — cos kx — cos fcj, 



(8) 



where fc^; and ky are the components of the wave- vector k. 
The mean field parameter Q and the chemical potential 
/i are determined by the self consistent equations 



S 



k 

Q = ^^ni^cosk^ 



(9) 



where Uk = [exp(wk/2^) ~ is the Bose occupation 
factor. Note that here the chemical potential /i is nec- 
essarily negative. At low temperatures, since the above 
sums are dominated by the modes Wk ^ T, Q ~ 5 and 

^ as T — > to ensure that the total particle num- 
ber remains fixed. A more precise analysis of ([9]) gives 
H ~ —T exp(—47r JS'^/T). The transition to the ordered 
state is explained by a pseudo Bose condensation of the 
bosons at the critical temperature^. 

Using the above results, we find the decoherence in the 
weak coupling regime is given by ^ with 



in 



^|A(k-q)p,5(£; + CJk-^^q) 



k,q 

X (rik + riq + 2n]f_nq) 



(10) 



where A(k) Aj exp(ik-i). In accordance with Refill, 

a multiplicative factor 2/3 has been taken into account so 
as to ensure that the spin correlation function satisfies the 
correct sum rules. In the zero temperature limit, since 
vanishes for uJk 3> T, the last term of PH]) is a Dirac 
function at i? = whereas the first two terms converge to 
a continuous function of E with the characteristic energy 
JS and a finite value, A(0)^/3J, at E = 0. For times 
t ^ (JS)^^, the sine function in ([7]) can be expanded 
leading to the usual short-time quadratic decoherence. 
For {JS)~^ < t < the last term of is essen- 

tially a Dirac function in ^ as it practically vanishes 
for \E\ >T. Moreover, the contribution of the first two 
terms of ((TU]) to Inr(i), -iA(0)2/3J, is negligible and 
then 



Inr(t) 



(11) 



where Aj^ = A(0) = X^i'^i- To obtain the decoherence 
for longer times, we evaluate T2{E) for energies \E\ <C T 

I 



using the approximations ojk — JS{k'^ 
q) ~ A(0). We find 



A 2 Tn2 

67rJ2S'2|^| 
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kl) and A(k 



(12) 



For times T^^ < i < lMr\ (US]) results in a smaU 
positive correction to the Gaussian decay (fTTjl . Inr(t) -I- 
4:{Af^St)^/3 ~ 2T(A/„jt)2 ln(Ti)/37rJ. For longer times, 
t ^ \fJ'\~^, the decoherence is Markovian with the rate 
r2(0). The non-Mar kovian corrections arising from the 
low energy behaviour of (|12p are logarithmic, Inr(t) -|- 
r2(0)i ~ (A/„T/7rJS'^)2ln(|^|t)/6. As anticipated, the 
rate r2(0) diverges in the zero temperature limit, im- 
plying that the asymptotic decoherence is Gaussian at 
T — 0. Note that the cross-over from the Gaussian de- 
coherence given by pT|) to the Markovian decoherence is 
extremely long at low temperatures. 

For the antiferromagnetic bath, following the steps 
outlined in Refs. and [ij, we obtain the following 
mean field Hamiltonian for the bath (up to a constant) 

^MF = Lk ^k(ak"k + ^k^k) where ak and /3k are lin- 
ear combinations of the Fourier transforms of the original 
operators aj and bi, respectively. The magnon dispersion 
is now given by 



Wk 



^k 



(13) 



where Ak = 2JQ{coskx + cos ky) and the corresponding 
self consistency conditions are 



nk+2 



(14) 



In the paramagnetic phase, T > 0, there exists a gap 
A = (/i^ — 16J^Q^)^/^ in the magnon dispersion for k = 
and the Neel ordering wavevector k = (tt, tt). An anal- 
ysis of dm) shows that as T ^ 0, Q Qo - 5* -I- 0.08 
and the gap vanishes as A ~ T exp{—2Tr\J\ps/T) where 
\J\ps is the spin stiffness of the bath^. The dimensionless 
parameter ps depends only on S, ps — 0.176 for S* = 1/2 
and Ps ~ 5^ for 5 > 1. 

A direct calculation of the dynamic spin correlation 
yieldai^ 



3/V2 



r/.2 



k,q 



e=±l 



^k^q 



+ e 



27T,kriq + ^k 



1 - e 



S{\E\+euji,-uj^) (15) 
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As in the ferromagnetic case, in the zero temperature 
limit, the sum of the terms oc n^nq is a Dirac func- 
tion at -E = whereas the other terms converge to a 
continuous function of E with finite characteristic en- 
ergy, \J\Qo, and value at E — 0, 2A(7r, 7r)^ps/3| J|(3o- 
Here also, for temperatures T <^ \J\S, this continuous 
function contributes to the decoherence only in the very 
short time regime t ^ (| JlS*)"^. To evaluate the decoher- 
ence for longer times, we remark that, for T <C \ J\S and 
<gC I Jl^, the sums over k and q in are dominated 
by the vicinities of (k, q) = (0,7r) and (k, q) — (7r,0) 
where tt = (tt, tt) permitting us to expand wj, and cjq to 
quadratic order in k and tt— q or vice- versa in conjunction 
with the consistent approximations |A(k— q)| ~ A(7r) and 
fi^ — A]f^A^ ~ 8| JjQo- Using these, we find the decoher- 
ence is Markovian for times t ^ with the rate 

r2(o) - ,^2^ 



37rJ2Q2A 



(16) 



where A^/m = A(7r). For times {\J\Sy^ < t < 
taking into account the terms with e = —1 in psp we 
obtain 



enhanced decoherence is observed in the vicinity of the 
transition temperature only if the central spin couples 
to the critical mode of the bath. We now show that this 
critical enhancement of the decoherence is not contingent 
on the absence of internal dynamics for the central spin. 
Due to the rotational invariance of H, any intrinsic dy- 
namic for the central spin can be described by the total 
Hamiltonian H' = H + £az/'2. In this mentioned 
earlier, the Markovian asymptotic behavior of the central 
spin state is characterized by two time a^^i^^ , a relaxation 
time Ti — r2(e)~^ and a decoherence time T2 given by 
T^'^ = Tf V2 + r2(0)/2. At the critical point, for any 
value of e, T2 vanishes and the resulting decoherence is 
faster than exponential. On the other hand, the behav- 
ior of Ti depends on the value of e. For |e| ^ \ J\S, the 
above study of the function r2 shows that Tj~^ reaches a 
maximum at T ~ |e|. This maximum grows as e — s- and 
we recover the critical enhancement of the decoherence 
for £ — 0, i.e., Ti = T2 — > 0. We finally remark that, in 
the low temperature limit, since Q is essentially constant, 
the spin environments we consider are practically similar 
to baths of independent but conserved bosons. 



\nr{t) 2± - ( Aa/m^t 



(17) 



As in the ferromagnetic case, the rate r2(0) diverges in 
the zero temperature limit and the asymptotic decoher- 
ence is Gaussian at T = 0. We note that the decoherence 
is faster as S" — > 00 i.e., when the spins become classical. 
In this case, since the associated scales fi,A 0, the 
Markovian regime is not accessible and the decoherence 
is Gaussian with a characteristic time oc S~^. 

We observe that the decoherence is qualitatively sim- 
ilar for ferromagnetic and antiferromagnetic interactions 
in the bath but with the important difference that Inr(i) 
is proportional to A'j^ = A(0)2 or to A^^^^ = A(7r)2, 
respectively. This can be understood as follows : an 



In conclusion, the two models studied in this paper 
clearly illustrate the phenomenal impact of fluctuations 
near a critical point on the decoherence and that to mini- 
mize decoherence it is essential to avoid continuous phase 
transitions in the bath. We expect our conclusions to be 
generically valid for any bath manifesting a continuous 
phase transition provided the central spin couples to the 
relevant critical modes. A natural extension of our work 
would be to study the decoherence and relaxation in- 
duced by a higher-dimensional bath in the ordered phase 
below Tc 7^ so as to compare the weakly fluctuating 
limit r ^ and the strongly fluctuating limit T T~ . 
Given the richness of the decoherence studied in this pa- 
per, it would be interesting to explore the decoherence 
induced by baths in quantum critical regimes. 
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